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By means of a formal expression of the Cornwall-Jackiw-Tomboulis effective potential for quark 
propagator at finite temperature and finite quark chemical potential, we derive the real-time ther¬ 
mal Schwinger-Dyson equation for quark propagator in Landau gauge. Denote the inverse quark 
propagator by A{p^) '/(— B{p^), we argue that, when temperature T is less than the given infrared 
momentum cutoff p^, A{p^) = 1 is a feasible approximation and can be assumed in discussions of 
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I. INTRODUCTION 

The Schwinger-Dyson (SD) equation for quark propa¬ 
gator is one of strong tools to research dynamical chiral 
symmetry breaking. A great deal of work on it in the 
zero temperature case has been made and a full demon¬ 
stration of chiral symmetry breaking in zero tempera¬ 
ture (^antum Chromadynamics (QCD) has been given 

0 However, the Schwinger- 
Dyson approach of chiral symmetry at finite temperature 
has not been fully explored. Such research is of impor¬ 
tant signihcance for understanding of the chiral phase 
transition at high temperature and even the color su¬ 
perconducting phase transition at low temperature and 
high density in QCD na , since the SD approach is a non- 
perturbative one which will give more reliable results. 

In present work on the SD approach of chiral phase 
transition, most are based on the imaginary-time for¬ 
malism of thermal held theory na and very few were 
seen in the literature which are based on the real-time 
formalism. In fact, in the real-time formalism, one can 
directly use some results in the zero temperature QCD 
consistent with the renormalization group (RG) analy¬ 
sis, thermalize them by thermal transformation matrices 
and obtain the expected expressions at hnite tempera¬ 
ture and chemical potential. In particular, the research 
in zero temperature case has shown that in the Landau 
gauge, the calculations can be simplihed greatly It 
is expected that at a hnite temperature, this advantage 
will be maintained if the Landau gauge is still taken. 

The SD equation for quark propagator can be directly 
derived from the extreme value condition of correspond- 
ing Cornwall-Jackiw-Tomboulis (CJT) effective potential 
[ 13 . In fact, if being limited to derive the SD equation, 
one needs only a formal expression for the effective po¬ 


tential rather than its a hnal explicit one. Conversely, the 
discussions based the SD equation of the mass function in 
quark propagator will be quite useful to simplifying the 
explicit calculation of corresponding CJT ehective poten¬ 
tial for quark propagator and this is also one of purposes 
of our research in this paper. 

The paper is arranged as follows. In SectE] we will 
give a formal expression of CJT ehective potential for 
quark propagator at hnite temperature and hnite quark 
chemical potential in the real-time formalism of thermal 
held theory and from which in Sect lllll we will derive 
in Landau gauge the SD equation for quark self-energy 
in hnite temperature and hnite chemical potential case. 
In Sect llVI an important possible approximation for the 
quark propagator will be explored and its practical sig¬ 
nihcance for chiral phase transition problem in QCD will 
be indicated. 


II. THE THERMAL CJT EFFECTIVE 
POTENTIAL FOR QUARK PROPAGATOR 

It is well known that at temperature T = 0 and the 
quark chemical potential /r = 0, the global havor chiral 
symmetries in QCD will be spontaneously broken owing 
to the vacuum quark-antiquark condensates {ip'if’) ^ 0 
j this makes quarks acquire their dynami¬ 
cal masses. For research on chiral symmetry restoring 
which could occur at a hnite T and p, the thermal effec¬ 
tive potential for quark propagator is essential. We will 
derive a formal expression of the ehective potential hrst 
by writing the CJT ehective action r[G, G*] in QCD for 
the quark propagator G and its conjugate G* at hnite T 
and p by 

rT[G,G*] =ri[G,G*] + r2[G,G*], (i) 


where ri[G, G*] and r2[G, G*] are respectively the con- 
*The project supported by the National Natural Science Founda- tribution from one-loop vacuum diagram and two and 
tion of China. more loop vacuum diagrams without quark self-energy 
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correction, since G in Eq.(l) is the exact quark propaga- where 
tor. In the real-time formalism of thermal field theory, 

we have S{p) = i/{^ + ie), S*{p) = —i/{^ — ie), = (6) 


ri[G',G*] = -zTr(ln 

-fTr(S'y + fTrl, (2) 

where St and Gt are 2x2 thermal matrix propagators, 
the superscript ”11” represents the 11 component of the 
corresponding matrix and the Tr is in functional sense. 
Assuming translational invariance, then in momentum 
space, we have M 

Srip) = MpS{p)Mp, Grip) = MpG{p)Mp (3) 


with the thermal transformation matrix Mp defined by 


Mp = 


cos0„ — e^^/^sin0„ 


gdM/2 sin 01, 


cos 0J, 


I3=l/T 


sin 0p = 9{p^)n{jf - p) + 9{-p^)n{-p^ + p) 


i(pO -p) = l/[e^(p“-'^) + 1 ] 


and 


Sip) - ^ , Gip) 


+ 1] 

(4) 

Gip) 0 \ 

(5) 

0 G*ip) j 


and 


Gip) = i/[Aip'^) Bip'^)+ie], 

G*ip) = -i/[Aip^) ^-Bip^)-ze]. (7) 

Gip) and G*(p) are the complete propagators of quarks 
with dynamical mass, where A(p^) and Bip^) have been 
assumed to be real functions. Thus in the momentum 
space, 

ri[G,G*] = -niNfN,{{tr[\nSTip)G^\p)]^^) 

+(tr[S'T^GT(p)]^^) - (trl)} 

= -llVi[G,G*] (8) 

where Nf and Nc are respectively the number of flavor 
and color of the quarks, H is the volume of space-time, 
tr only represents the trace of spinor matrices and the 
denotation (• • •) means the integratioii fd'^p/(27r)‘^. By 
means of Eqs.(3)-(6) and the relation [l2| 

In St ip)G^\p) = MplnlSip)G-^ip)]M-\p), 

we obtain from Eq.(8) the one-loop effective potential 
Vi[G, G*] corresponding to ri[G, G*] 


Vi[G,G*] = lA^/A^c ((cos^0p{trln[S'(p)G ^(p)]-ftr[S' ^(p)G(p)]}) 

-f(sin^0p{trln[S'*(p)G*"^(p)] -f tr[S'*"i(p)G*(p)]}) - (trl)) . (9) 

In the real-time formalism of thermal field theory, the interaction Lagrangian between the quark fields '0 and the 
gluon fields A“(a = 1, • • •, 8) can be expressed by 


a=l r=l,2 

where r = 1 and r = 2 correspond to physical and (thermal) ghost fields respectively, A“(iV^) is the generator matrix 
of the color gauge group S'Gc(3) in the representation JV^ of the quarks and go is the bare gauge coupling constant. 
Assume that only the two-loop vacuum diagrams are included in r 2 [G, G*], then we can write in the momentum space 

r2[G,Gl = -^((5o"Tr7'^A“(iVJG“(p)7"A'’(iVJG“(g)[i^;r,*'(p-g)]“))Ar/fi 

+ ^((5oTr7'^A“(iVJG)?(p)7"A'(iVJGii(g)[iA;)^(p - <?)])?))Ar^ll 
= -nV2[G,G*] ( 10 ) 


where [D'^^ip — q)]T is the complete thermal gluon 
matrix propagator, ((• • •)) represents the integrations 


f f dApd'^q/i2'K)^ and Tr is now only for the quark- 
dependent matrices (flavor, color, and spinor etc.). After 
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the renormalization procedure is implemented, we can 
replace by the RG-invariant running gauge coupling 
\{p—gY \) and \D'^^{p—q)\rp by the tree diagram gluon 
propagatoiJZ?“(^(p — q)]rp, i.e. we will take the approxi¬ 
mation 

5o [D';^\p - 9)] T ^ 9 ' [(P - 9 )"] [DfAP - 9)] (11) 


with the thermal transformation matrix defined by 


Mk 


f cosh 0 fe sinh 0 ^ 
sinh 0 fc cosh 0 fc 


n(fc°) 


1 

6^1'=”! - 1 


sinh 0 fe = sjn{k°), 

(13) 


In Eq.(lO) the complete thermal vertices (p,q){r = 
1 , 2 ) have been replaced by the tree diagram vertex 7 ^^. 
This is because when A(p'^) and B(p^) in the inverse 
quark propagator G~^(p) are real functions, the Ward- 
Takahashi (WT) identity at finite T and p (if the ghost 
effect in the fermion sector is neglected) essentially iden¬ 
tical to the one at T = ^ = 0, thus in the Landau gauge 
we can make the couplings between the complete vertices 
q)(r = 1,2) which submit to the WT identity 
and the thermal gluon propagator are equivalent to the 
one between the tree diagram vertex 7 ^ and the thermal 
gluon propagator 0 The tree diagram thermal gluon 
matrix propagator [D^j^{k)]T can be expressed by jl^ 

[DtM]^ = [D^..{k)]^ = MkD^,{k)Mk 

( 12 ) 


and 

= ( "’"f > oyt)) ■ (»» 

In Landau gauge. 


Df,„{k) _ Ti f kf,k^ \ 

D%{k) k^±ie)- 


(15) 


By means of Eqs.(3)-(5) and Eqs. (11)-(14), we may 
obtain from Eq.(lO) the two-loop effective potential 
V 2 [G, G*] corresponding to r 2 [G, G*], expressed by 


E 2 [G,G*] = ^-NfN,C2{NX{9%p-gfMrG}r^{p)YG}r\q)\ [D^,{p - q)]]^)) 

-^7V/WG2(iVj((52[(p_ g)2]tr[7^G^2(p)7'^G^i(g)] [D^.{p-q)]]:^)). (16) 

with 

G^{p) = cos^ 9pG{p) — sin^ 6pG*{p), G^{p) = — cosOp sinOp e^^^‘^[G{p) + G*(p)] = —e^'^Gy (p), (17) 

[Dp„{p-q)]l^ = Dp„{p-q)-gf,^27Tn{p°-q°)S[{p-q)‘^], [Dp„{p-q)]l^ = -gp„e/^^P'‘-'^°^^‘^2Trn{p°-q°)6[{p-qf] (18) 

and 


lVcG 2 (fVJ =^tr[A“(fVjA''(iVJ], G 2 (iVj 

a,b 


2N, 


The total effective potential will be 


E[G, G*] = El [G, G*] -f E 2 [G, G*] 


(19) 


III. THE THERMAL SCHWINGER-DYSON 
EQUATION FOR QUARK SELF-ENERGY 


The effective potential E[G, G*] is in fact the func¬ 
tional of the thermal propagator [G{p)]rp and 
q)]y (rs = 11,12). The Schwinger-Dyson equation for 


quark self-energy will come from the extremum condition 


SV[G,G*] 
6G^ (p) 


( 20 ) 


In view of Eq.(17), equation (20) is equivalent to 


6 

6G{^ 


S 

5G*{p) 


U[G,G*] =0 
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which, by using Eqs.(9), (16) and (19), can be reduced 
to 

E(p) = [1-A{p^)]^ + B{p^) 

= -^C2m{g\p-qf]^^G}rHqh''[D,Ap-q)]¥) 

( 21 ) 

Equation (21) is just the SD equation of the complete 
self-energy E(p) for single flavor and single color quark 

I 


to one-loop order at finite T and p. It is noted that the 
second term in the right-handed side of Eq.(16) has no 
contribution to Eq.(21). By acting and tr on the two 
sides of Eq.(21) separately, we can obtain the equations 
of A{p'^) and B{p^), 


A{p^) = 1 + *^£l(52[(p- ,)^]tr[^7'^G“(<z)7l [D,^{p - q)]]^) 


( 22 ) 


B{p^) = -i [(p - qf] tr {q)Y] [Df,u{p - q)]^^) ■ (23) 

Erom Eqs.(7) and (17) we have 

Gt (<?) = [A{q^)^ + B{q^)] -g ~ 27rsin^ egS[A^{q^)q^ - B^{q^)] | . (24) 

By means of Eqs.(18) and (24), and the result that 

tr(|l 7 ^ 4Y)D^^{p- q) = -AE{p,q) 


with 


Eip, q) = l- - 


- ■ - ■ > /(P-,)t 


p^ + q^ + 


[p - qf 


(25) 


equation (22) can be changed into 


AiP) = l + 


^ 2 (q, 2 ^q 2 _ _|_ jg. 

x{ - iE(j), q) + ATr{p ■ q)n{p° - q°)S[{p - qf]'^y 
It is obtained that the integration in Eq.(26) 

Eip, q) 


27rsin^ 0q6[A'^{q'^)q'^ - B'^{q'^)] | 


(26) 


g^[ip-qr]A{q^) 


A^{q^)q^ - B^{q^) + ie 


d-^q 


g%p-qr]Aif) 


E{p, q) ^ „ 

A\f)f+B\f) ’ 


where we have made the changes of variables = ip^, p* = P* (* = 1, 2, 3) and the same to q ^ q after the Wick 
rotation, and the assumption that 

g^ [(p- - qf] = 0(p-2 - f)g\p^) + e{f - p^)gHf), 

and then consider the fact that 

J dnqE{p, q) = 0. 

As a result, equation (26) is reduced to 


A(p2) = l-27r^^^^(^g^[{p-qf]Aiq^){Eip,q)smHqS[A^iq^)q^-B^iq^)] 

+ ( y^ 2 (g 2 )g 2 _^g 2 (g 2 ) + ^g + *27rsin^ eq6[A^{q^)q^ - B^q^)y ■ 2(p • q)n{p° - q°)6[ip - ?)"] })• (27) 
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On the other hand, by using Eqs. (18) and (24), equation (23) can be changed into 

X{ (p _"jt + fe - 8""(P° - g’ldfr -?)*]})- (28) 

Equations (27) and (28) are the required coupled equations of the function AfjP') and B{p^) in the quark self-energy 
S(p) in Landau gauge in the real-time thermal QCD. It is indicated that when T = /i = 0, all the thermal correction 
terms can be removed, hence equations (27) and (28) will be reduced to 

= U B(P-) = -aC2(iU){g^l(p - AU ' (P- <!y + ie ) 

which are just the SD equations for quark self-energy functi on in Landau gauge at T = ^ = 0. _ 


IV. DISCUSSION OF A{p'^) 

Equation (27) shows that, at a finite T and A{p‘^) 
has extra thermal correction terms to unity. Based on 
the consideration of dimension, A(p^) can be generally 
written by AijP') = I — (T'^/p^)F, where F depends on 
only the dimensionless variants consisting of T, p, 
etc.. Assuming F is finite, then it can be deduced that 
at a high and a low p^, the thermal correction terms 
would be not small and A{p^) = I would not be a good 


approximation. However, we will argue that for some 
physically interesting range of T, the thermal corrections 
of A{p^) are negligible hence the approximation A(p^) = 
1 may still be valid. 

Denote the quark mass function by 

m\q^)=B\q^)/A\q^)2 (30) 

then A(p^) can be expressed by 


A(p 2 ) = 1 - 


+5' 


C2(iVJ f dU 1 


'<»>(; 


(27r)3 ^(< 7 ^) 
1 


<lf]Eip, q) sin^ 0gS[q'^ - m'^{q'^)] 

- z27rsin2 - m2(g2)]^(p2 q-^)n{p° - q°)S[{p - qf] }. (31) 


+ ie 

Let q'^ = ml is a real root of the equation = m^(g^), i.e. mi obeys the equation 

2 2 / 2 \ 
mi = m (mi), 


(32) 


then we will have 


6[q^-m^iq^)] 


/K) 


Sy^-mj)], f{ml) 



dm^{q^) 

dq^ 



(33) 


In fact, one can identify m{q^) with the running quark mass function consistent with the renormalization group (RG) 
in QCD 13, then it will be a descent function of q^ and we always have f{m\) > 1. To estimate the order of 
magnitude of the thermal corrections of A(p^), we will take the approximation 


gy{p-qf]E{p2q) e{\p^\-y]) g^p^)- 




2 \ 2 -| 


+0(i9^i-ip2|)y(a 1 - 


(p 2 - q^) 


2 \ 2 -| 


1 1^1 1 
A{q^) q^ — nn?{q'^) + ie A(mi)/(mf) q^ — m\+ie 


and 
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^T^T'^h{yi,r) 


with 


and 


we can write 


where 


1 /•“ 

-f 3 (j/i,r) =dx 


0 \/x‘^ +yl 


o^x‘^+y\-r _|_ 


+ (-r ^ r) 


mi y 

yi = r= - 


with 


lip^) = 


1 


/(to?) 

3 ^( 0 ) + m\ 


^f{p^)h{yur)+ g 


J d^qn{q^) 6 {q^) = ^TT^T^, 

Aip^) = l-a{p^)/Aiml), 

g^ioy 


(IpP +p 2 )£(p 2 ) 

p"^ -ml + 2xT(jP + 1^) 0 


T 


1 - '»! -r yiJ -r |jv|; , n _ qn 

127r2 (|p|2 f,x_iY p"^-m\^-2xT{p^-\p\) J 


2 l„2|^ 2/ 2 n P 


~g\py = e{mi-\py)g\mi) 


2{ml+pl) \ m 


mr / ^ ml 




2 p 2 


(34) 


(35) 


(36) 


In Eq.(35) we have introduced the infrared momentum 
cutoff pI by the replacement 1/p^ ^ l/(bP +4'?)^(p^) 
with e(p^) = and in Eq.(36) by the replacement 

l/2mi ^ l/2{m'l+pl). Whenp^ = mf, both 5 ^(p^) and 
the integral with the logarithm functions in Eq.(35) are 
equal to zeroes, hence we get 


a(m?) 


gHo)T^ 

9/(mf)(m? +p2) ■ 


(37) 


When 7 ^ ml and if |p^ — ml\ > is assumed, we 
can expand the logarithm functions in Eq.(35) in power 


of 2xT{p^ ± |p|)/(p^ — ml) (noting that the main con¬ 
tribution to the integral over x comes from the region 
of a; < 1), then in the approximation with the infrared 
momentum cutoff pi included that 


(p 2 — rnlY 


= 0{ml - Ip^I) 


-{ml+pl)\ 


+B{\p^ \ - ml) 


{p^y 


we can obtain that 


a(p^) = e{ml-\p^\)2a{ml) 


-|-6>(|p^| - m?)2a(m?) 


1 - 


\py J Sg^ml) 

|p2| +pI\ 27r2g2(0) 

TO? +pI I 3g^(p^) 

' -hp2)£:(p2) I 27r2p2(0) 


h{yi,r) - 1 - 


47r^ 3p° 


-f 


T2 


15 


1 -^ /3(.r,r)I^ + l + 


ml + pi J 

47r^ 3p°^ +p^ T'^ 


15 


pZ 


pZ 


(38) 


From Eq.(34) we can get the equation obeyed by A(ml) which has the solution 


z4(m?) = 1 — a(m?)/z4(m?) 


2 i(TO?) = -{l + [l- 4a(m?)]i/2}. 


(39) 
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When temperature T is low enough so that 4a(TO^) < 1, 
A{mi) will have the order of magnitude of one, and it 
is seen from Eq.(38) that we always have a(p^) <C 1, 
whether \p^ \ <C rn\ or \p^\ ^ m\. Hence we can obtain 


H(/) 1, 


9/(m? 

1 +pI ^ V(o) ■ 


(40) 


In this case, A{mi) becomes complex and we may write 
the module of a{p'^)/A{ml) by 


ajP^) ^ Hp^)\ _ Hp^)\ 

A{mf) ll+iy'Aaiml) - l|/2 


When temperature T is neither very low nor very high 
so that 4a(mf) > 1 but T'^/{m\ + p^) < 1, the series 
expansion in T'^/{m1 + pi) in Eq.(38) may still be used. 


Noting that the terms contained in the brace in Eq.(38) 
have the order of magnitude of unity, the module may be 
approximately estimated by 


Hp^)\ 

jl/2(.^2) 


— \p^\)2a}/'^{m\] 


|p2 


m+pl 


■ e{\p‘^\ - ml)2a^/^{ml) 


mt 


i\p\^+piy 


(42) 


In spite of 2^ya{mf) > I, when |p^| > mf, 

|a(p^)|/a^/^(mf) <C 1 is obvious if p^ is large enough. 
On the other hand, when |p^| < m\, we may make a nu¬ 
merical estimation of \a{py \/in practical case 
of QCD. If the running gauge coupling g^{p^) in QCD 
with the infrared momentum cutoff pi is taken as 




'■QCD 


1 is not very large and it is easy to verify that it will 
decrease as T goes down and the infrared momentum 
cutoff p1 goes up. It is also seen that the value of m\ 
either lower or slightly higher than the infrared momen¬ 
tum cutoff p1, so the region in which |p^| < m\ is quite 
limited, and it may be presumed that the deviation of 
A(p^) to 1 will cause only quite small effect to the total 
results. Hence we can neglect such deviation and take 
A{p^) ~ 1 for all the value of p^ in the conditions 


where A = 24/ (33 — 2Nj) and Aqcd is the RG-invariant 
mass scale parameter, then the quark mass function 
m{p^) (when only dynamical quark mass is involved) can 
be expressed by 


9/(mj) 

4g2(0) rril^-pl 
To sum up, we can assume that 


(43) 


V 3 J \p^ +p1J \ A2 ] 


m 


Irl+P: 


2 , A-2 


A2 

^^QCD 


where (j) is the RG-invariant quark-antiquark condensates 
[i|. Thus Eqs.(32) and (33) lead to 


mi = 


27r2A \(j)\ 


m\ + p1\^A-^ 


^ ^^QCD ^ 


and 


f{ml) = 1-b 


2m\ 
m? -I- E 


1-b 


2-A 


(>“M 

^ ^^QCD ^ 


2n -1 


It may be checked that, taking p^ = mi/2, when Nf = 3 
and tc = \n{p1/Al^cD) = O.I, for |/>|/A|c£, = 0.1, we will 
have m\lA qcd — 0.64 and obtain |a(p^)|/a^/^(mi) = 
0.5 and 0.86 respectively for T/{m\ pi) = 0.5 and 
T/pc = 1; for 1(1)1/A^cd = 0-5, we will have mi/A qcd ^ 
1.2 and obtain |a(p^)|/a^/^(m2) = 1.33 and 1.75 respec¬ 
tively for T /(mf -I-Pc) =0.5 and T/pc = 1. These results 
show that when |p2| < jn\, A(p2) ~ 1 is not a good ap¬ 
proximation. However, the difference between A{p^) and 


A(p^) ~ 1, for 0 < < m-i -I-Pc (44) 

is a feasible approximation. When only dynamical quark 
mass is involved, the minimal mi will be zero, and when a 
current quark mass is also included, the minimal mi may 
still be quite small, so the allowed highest temperature 
limited by Eq.(44) may be reasonably supposed to be 
T = Pc = 6 *“/^Aqcd = l.OhAQcn (e.g. for = 0.1) 
and will increases as a larger pc is taken. Hence, if the 
chiral phase transition temperature is lower than Aqcd 
\\‘A Il9| , then the above arguments indicate that within 
the whole range of temperature related to chiral phase 
transition, A(p^) = 1 is a feasible approximation and 
may be used in the discussion of chiral phase transition 
problem in QGD based on Schwinger-Dyson approach. 
This will certainly greatly simplify calculations in this 
nonperturbative method. 
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